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ABSTRACT 


Rapid assessment of littoral oeean physieal and ehemical conditions has 
represented a great challenge in recent years. For the physical part, time constraint may 
limit to using the barotropic mode. But because rain can be significant in St Andrew Bay 
system, Florida, fresh water, even when rivers lack, is a prevailing salinity regulator 
through ground seepage. Therefore, studying the baroclinicity and forcing mechanisms 
should provide a reasonable guidance help for deciding which will be ignored or included 
when modeling the area of interest. For the chemical part, considering either national 
security relevance or shipping hazard, a release of chemicals used in nearby offshore oil 
platform stands in the domain of possibilities. 

A coupled hydrodynamic and chemical model, WQMAP and CHEMMAP 
developed at the Applied Sciences Associates, is used in this study with WQMAP for 
evaluating baroclinicity and forcing mechanism impacts and with CHEMAP for 
estimating the consequences of a hypothetical release of chemical. A stochastic model 
was applied for determining probable distribution and concentration resulting from the 
release at a given location. 
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I. 


INTRODUCTION 


The U.S. Navy has requirements for the prediction of currents and tides 
worldwide in geographical basins ranging from open ocean to continental shelf, estuary, 
and rivers. The Naval Oceanographic Office (NAVOCEANO) implements and runs 
several types of 2D and 3D circulation models for many geographical domains around 
the world. Most of these are run in operational modes and generate daily products for the 
Fleet. Although the mission is to provide nowcasts and forecasts worldwide (deep and 
shallow water), many operations take place in coastal bays and rivers. The model 
products in these domains support a wide variety of applications ranging from Special 
Operations, Mine Warfare, Expeditionary Warfare, Object Drift, Search and Rescue, 
Chemical and Oil Spill. St Andrew Bay system, hosting the Naval Surface Warfare 
Center, is a perfect theatre of this study because it has the intra-coastal highways, and is 
in the neighborhood of offshore oil drilling platforms. Rapid environmental assessment is 
important for this region. 

High resolution, 3D baroclinic nested littoral-open ocean models are usually used 
to assess complex littoral zone oceanography with forcing functions. The purpose of this 
study is to determine if a 2D barotropic, un-nested littoral circulation model can satisfy 
certain Navy applications in coastal bays dominated by tides and local wind forcing 
within time scales ranging from “instantaneous” (what a diver will feel in the water) to 
days (the drift of floating objects). Savings in model implement, tuning, and validation 
are of outmost importance when using rapid response systems as authorities requires 
shrinking time delay between accident occurrence and the first decision. 

Applied Science Associate, Inc. (ASA) out of Narragansett, Rhode Island has 
developed a three-dimensional, time dependent hydrodynamic model, the Water Quality 
Mapping and Analysis Package (WQMAP) coupled with rapid response model for oil 
spill tracking (OILMAP), search and rescue operations (SARMAP) or chemical release 
(CHEMMAP). The possibility of running these hydrodynamic models either in 2D- or 
3D-version, the fact they have been used for years by NAVOCEANO, increased the 
significance of their use. 
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The lack of in-situ measured data did not constrain this study as we focused on 
the effects each forcing mechanism induced to the system. We hence built up a control 
run as a reference picture for circulation pattern, current speed and vertical salinity 
distributions. A careful scrutiny of local forcing pinpointed their ability in influencing the 
circulation in St Andrew Bay system and the consequences of each forcing mechanism 
on the system were examined. We finally choose to run the hydrochemical model 
CHEMMAP as rapid response system to evaluate the response of St Andrew Bay to a 
hazardous chemical release. 
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II. 


ST ANDREW BAY OCEANOGRAPHY 


The continental shelf of the Gulf of Mexico includes the broad Western Florida 
Shelf that extends to 200 km in width (Roberts et ah, 1999). Major relief along the shelf 
in this region is associated with roughly shore-parallel ridges, related to the positions of 
former shorelines during the rise of sea level after the latest Pleistocene glacial 
maximum. The northeastern Gulf of Mexico continental shelf off the panhandle of 
Florida, Alabama, and Mississippi ranges in width from 25 to 125 km, and depths at the 
shelf break range from 60 to 100 m. 

Shelf circulation in the Gulf of Mexico exhibits distinct patterns in different 
regions. Off western Florida, the circulation is complex, influenced largely by wind and 
Loop Current forcing. In the northeastern Gulf of Mexico, flow appears to be responsive 
to wind forcing. Average flow along the inner shelf is onshore near the bottom and 
offshore in the upper water column (Wiseman and Sturges, 1999). 



Figure 1, St Andrew Bay system 

(http://science,gulfcoast,edu/lindafitzhugh/St,%20Andrew%20Bay,jpg) the 20-May- 

05 
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A, ST ANDREW BAY SYSTEM 

The St Andrew Bay system is loeated in the western part of Florida State elose to 
Pensaeola. It is a long, narrow and rather shallow (water eolumn extends from 4 to 10 m) 
bay oriented in the northwestern direetion. Deer Point Lake loeated at the northern edge 
of North Bay, provides the major freshwater inflow into the estuary, along with a number 
of smaller ereeks. Two major passes. East Pass and West Pass, have provided surfaee 
water eonneetions via its eentral part (St Andrew Bay) with the Gulf of Mexico. West 
Pass was artificially cut in 1934 in Shell Island as the primary navigation channel to the 
Gulf, while most exchange between the estuary and the Gulf had historically occurred 
through East Pass. East Pass was recently filled in (1998), however, by the movement of 
shoreline sediments. Dredging was conducted for reopening the pass in 2001. None of 
these channel widths exceed 1 km. St Andrew Bay system also links the Intracoastal 
Waterway on both the very end of the East Bay and the West Bay. The waterway depth is 
constantly maintained at 4m which causes some water exchanges to occur at these 
locations. 

Drainages into these bays are small and mainly due to groundwater seepage. Eor 
example, it is suggested that almost two thirds of the creeks which flow into Deer Point 
Take are supposed to be spring fed (Musgrove et ah, 1965). Nevertheless their influence 
can be crucial when modeling the whole system. 

In terms of estuarine classification, the four basins are generally positive, i.e., 
drainage inflow exceeds evaporation (Pritchard, 1952) except in St Andrew bay where 
neutral conditions, due to its connection with the Gulf, are found (Ichyie and Jones, 
1961). The bottom contours are part of the WQMAP geographic information system 
(GIS) layer and have been provided by NOAA. They depict a 5.2m average depth for St 
Andrew bay, the mean depth respectively being 2.1, 1.7 and 2m for East, North and West 
bays. 
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Figure 2, St Andrew Bay bathymetry 


Temperature is rather uniform within the water column through out the bay 
system with very small vertical variations even in summer (the depth is too shallow for 
the thermocline to really develop). But its seasonal variability is evident. A five year- 
study of surface temperature from both inside and outside the bay, conducted using 
NOAA time series data, shows that the mean temperature varies between 15“ C in winter 
and 30“ C in summer. The comparison of these two 5-year-long time series (between 
2000 and 2004) portrays a strong zonal sensibility and the difference ranges from 0 
during winter time up to 2.5“ C during summer (Fig. 3). The rather shallow bay allows 
the water column to heat up very quickly; the maximum difference is reached no later 
than March and lasts almost 9 months. 
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Figure 3, Temperature difference between St Andrew Bay and Gulf inner shelf 

Salinity, however, presents some zonal features. Close to the Gulf of Mexico 
connection, the salinity is almost constant within the water column (gradient of 0.17 
psu/m) and reaches its highest value (34.3 psu). Its value can drop down to 9 psu at 
surface nearby with a gradient of 3.3 psu/m (Ichyie, 1961; Blumberg and Kim, 2000). 
Nethertheless, the 24 h variation of salinity does not correlate closely with variations of 
the tidal current, specifically away from the bay entrance. The salinity follows the net 
flows caused by wind driven currents and freshwater inflow events. The hydrographic 
study conducted by Ichyie (1961) proved the water column stability being time 
independent and mainly subjected to salinity vertical variations (Sverdrup et ah, 1942). 
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A, TIDES 

The Panama City Beach (30° 12.8’N - 85° 52.7’W) and the St Andrew Bay (30° 
OO.l’N - 85° 40.0’W) data study shows that according to the form ratio (ratio between 
diurnal -Kl+01- and semidiurnal -M2+S2-), the tidal constituents are diurnal (ratio 
greater than 3). High and low water times study show the tidal waves first enter the bay 
through East Pass (5 minutes difference between East and West Pass) which correlates 
the general motion of tidal waves. The maximum current is found in the West Pass and 
can exceed 0.5m/s. It then decreases as low as 0.15m/s at the western end of West Bay. 
Most of the time, ebb currents duration is longer than flood current duration, mainly due 
to the effect of drainage, source of a quasi permanent outward current in the upper layer 
if any (Pritchard, 1952). This process causes “tidal pumping” phenomenon (Eischer et al, 
1979). 

Compound tides, often important in shallow water and hence on inner shelves, are 
generated by the nonlinear interactions of the primary constituents. The long-period tidal 
constituents are not resolvable in the month-long series, and contributes little (1.1 cm) to 
the tidal elevation. It is crucial to understand any seasonal effects which may occur on 
tides. Tidal results from a yearly sea-level time series were compared with those from a 
monthly subset of the sea-level time series to examine the seasonal effects. Additionally, 
this provides insight into the effects of unresolved tidal constituents contained in shorter 
time series. Overall, none of the major diurnal and semidiurnal tidal constituents differed 
more than 1 cm between the year- and month-long series except K1 constituent which 
varies significantly. 

The five main constituents have been computed though 12 successive month-long 
time series collected in 2004 outside St Andrew Bay and are associated with a 95% 
confidence interval. This general trend is consistent with the observations of the last 5 
years and denotes the occurrence of a seasonal dependence of both K1 and S2 
constituents. A study, conducted in 2003 by Alvarez (personal publication), shows that 
sea breezes may be considered as a major source of the seasonal variability in the K1 
tidal dynamics in near-coastal shallow regions this fact being explained as follows: if the 
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K1 tidal response to boundary foreing is subjeeted to the sea-breeze impaet, the 
superposition of the K1 signal and the SI signal, independent of what is its origin, will 
produee a seasonal modulation of the K1 tidal dynamies and, henee, seasonal variations 
in the K1 tidal eharaeteristies. This finding results from the well-known expression for an 
amplitude- and phase-modulated oseillation. This statement is supported by experiments 
eondueted in the Bay of Cadiz. However, the SI eoeffieient does not show up during a 
monthly eomputation (the frequeney step is too small for separating SI from Kl) and 
only portrays an amplitude of 0.3om (2% of the Kl eonstituent) in the year-long time 
series study. It is also to notiee that the eonfidenee interval is very large during 
September as the data were impaeted by a series of storm and a hurrieane. 



2 4 6 8 10 12 

Months 


Figure 4, Diurnal constituents seasonal variability (at Panama City 

Beach) 
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Figure 5, Semi-diurnal constituents seasonal variability (at Panama City 

Beach) 

The only difference between both tidal analyses is the existence of residual 
semidiurnal energy, which accounts for small, erratic oscillations in the residual of the 
month-long series. While these residual oscillations are small compared to the dominant 
tidal constituents, it illustrates the problem of tidal analysis resolution inherent in short 
record lengths. 

B, FRESHWATER SUPPLY 

Freshwater discharge in the Gulf of Mexico totals approximately 1110 km per 
year and is dominated by the Mississippi and Atchafalaya rivers, which contribute 55% 
of this discharge (Solis and Powell, 1999). The extent of the impact of this discharge, 
both vertically and horizontally, is variable and modulated by discharge rate and 
processes that disperse the plume, primarily wind stress. The river water accumulates in 
the inner shelf, forming a low salinity band that supports a persistent baroclinic flow to 
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the south (Atkinson et ah, 1983). The mid shelf zone (21-40 m isobaths) is vertieally 
homogeneous during the fall (September, October, November) and winter (December, 
January, February), due to enhanced wind mixing and decreased runoff. This region 
undergoes a transition to a vertically stratified state in the spring (March, April, May) and 
summer (June, July, August) as wind mixing decreases and runoff increases (Atkinson et 
ah, 1983). In addition to runoff, groundwater sources have been shown to be important 
in the SE region of United States and prevail in the St Andrew bay system. 

C. WAVE CONDITIONS 

Waves largely correlating the wind direction, they primarily come from the east 
(23.5%). 2004 measurements from buoy 42039, located at 28.80N-86.06W, largely 
confirm the average computed from 1996 to 2001, the correlation coefficient being 0.77. 



Direction [degree] 

Figure 6, Wave direction comparison between 2004 and 1996-2001 average 
(buoy 42039, located at 28.80N-86.06W) 
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The monthly significant wave height (SWH) and wave peak period averages in 
figures 7 and 8 enhance the rather calm period between April and August. The yearly 
SWH mean is 1.1m but its median stands at 0.86m showing clearly the impact of 
hurricane and tropical storm period. The comparison between data collected in 2004 and 
1996-2001 period averaged data do not depict any difference and present a 0.998 
correlation coefficient. 


42039 SIGNIFICANT WAVE HEIGHT ( METERS) 12/1995 - 12/2001 



Figure 7, Average wave height between 1996 and 2001 
(http://www,ndhc,noaa,gov/images/cHmplot/42039_wh,jpg) the 20-May-05 
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42039 DOMINANT WAVE PERIOD { SECONDS) 12/1995 - 12/2001 



Figure 8, Dominant wave period between 1996 and 2001 
(http://www,ndbc,noaa,gov/images/climplot/42039_wp,jpg) tbe 20-May-05 



Figure 9, Significant wave beigbt comparison between 1996-2001 and 2004 
data sets (buoy 42039, located at 28,80N-86,06W) 
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However, wave generation in St. Andrew Bay is limited by fetch. Local harbor 
authorities state that wave height at Dyers Point is approximately 1.5 to 2.1 m during 
southeasterly hurricane force (>33m/s) winds. Winds from other directions would result 
in lower wave heights. Wave heights may reach 6 m in the outer portion of the entrance 
channel, but do not reach the inner part of St. Andrew Bay. 

D, ALONGSHORE CURRENT AND LITTORAL TRANSPORT 

The main current affecting the surf zone is the alongshore current created by 
waves breaking at an angle to the shore. The magnitude of the alongshore current 
depends on the breaking wave characteristics, breaking angle and local bottom and shore 
configurations. 

The alongshore currents are responsible for sand transport along the coast. For the 
study area, the net littoral transport is generally westward as the predominant waves are 
from the east and is estimated to be about 65,000 cubic yards per year. This estimate 
seems to be in agreement with the field evidence as indicated by the lack of strong 
erosion and accretion at the west and east side of the gulf entrance to the St. Andrew Bay 
(U. S. Army Corps of Engineers, 1971). However, the accumulation of hurricane-induced 
waves during the 1990’s is likely to be responsible for East Pass closure that happened in 
1998. Van de Kreeke (1990) showed that a two-inlet system could not be stable for a one- 
bay system. Even without erosion evidence, one of the inlets has to close. Moreover, Jain 
and al (2004) showed that both inlets were actually unstable and hence will request 
continuous dredging maintenance unless St Andrew bay is enlarged. 

Bay County, to which Panama City belongs, has been scouting the shoreline since 
1971 when Elorida State made each county responsible for providing data on erosion 
trends, vegetation line, bathymetry surveys. Prom 30 years of data, the erosion causes 
have been split into two main processes. 

1, Long-term changes 

All beaches are experiencing erosion, the more severe locating on the northern 
edge of main entrance channel with a rate of 2.44 m/yr. This process is associated with 
the construction of St Andrews Inlet in 1934. Similar rates occur along Shell Island and 
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vary from 0.53m/yr in its western part to 1.95 m/yr elose to the secondary channel 
entrance. The lack of construction on Shell Island does not make the eroding process as 
critical. However, bypass placements of sand have been helping in the stabilization of the 
shoreline. 

2. Short term changes 

These changes are mainly due to extreme weather system hitting the area. 
Hurricane Eloise in September 1975 caused a 100,000 m beach-dune erosion. Hurricane 
Opal in October 1995 impacted the area by dredging almost 100,000 m from both 
beaches and dunes (Bureau of Beaches and Coastal Systems Division). Finally, during 
hurricane Ivan event in September 2004, beaches alongside Bay County suffered a 
vertical loss of sand between 1.2 and 1.5 m (Bureau of Beaches and Coastal Systems 
Division, 2004). 
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III. 


ATMOSPHERIC FORCING 


A. WINDS 

Located at average latitude of 30“10’N, the St Andrew Bay system is under the 
influenee of rather weak winds. The elimatology of wind foreing in the Gulf of Mexieo 
exhibits a pattern of seasonal variation. During the fall and winter months (mid- 
September to mid-February), winds are primarily from the north. In the late spring and 
summer (from Mareh to July), the northern Gulf of Mexieo is dominated by tropieal 
weather with winds mainly southerly. During summer, the influenee of the subtropieal 
high (Bermuda High) inereases as the frontal zone between subtropieal and mid-latitude 
air masses moves north and out of the Gulf. Weaker pressure gradients and, henee, 
ealmer winds assoeiated with high pressure produee less vigorous wind stress foreing of 
oeeanie or, in partieular, shelf eireulation. During summer, warm fronts move generally 
from south to north. The average data eolleeted throughout the last 60 years show very 
uniform winds. 



Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

Dir. 

N 

N 

SSE 

S 

S 

WSW 

WSW 

E 

ENE 

N 

N 

N 

Speed 

3.5 

3.5 

4 

3.5 

3 

3 

3 

2.5 

3 

3 

3 

3.5 


Table 1, Average wind direction and speed (m/s) 


B. RAINFALL 

The preeipitation reeords from the last 30 years at Panama eity airport eorrelates 
the maximum rainfall with the northern (winter) and southern (summer) winds. However, 
beeause of their tropieal influenee, rainfalls are often made from showers whieh may vary 
between two loeations inside the bay, this phenomenon being pinpointed in former study 
(Blumberg and Kim, 2000). 



Jan. 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

rain 

2.18 

1.89 

2.42 

1.46 

1.48 

2.51 

3.35 

2.90 

2.52 

1.39 

1.76 

1.53 


Table 2, Average rainfall (cm) at Panama City Airport 
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C. STORMS 

The low wind speed average values must not hide episodie storm events this 
region has to faee. The outstanding feature of the U.S. Gulf Coast region is its loeation on 
the north shore of the Gulf of Mexico and its orientation perpendicular to normal tropical 
cyclone tracks as they move more or less northward out of the tropics. Also of 
importance is the region's position between 25 and 30 degrees north latitude. This is 
within the normal location of tropical cyclone curvature change, which oscillates between 
latitudes 25N and 35N during the tropical cyclone season. This latter factor is significant 
since it is the character of tropical cyclones to slow and intensify during this stage. As an 
example, last September, hurricane Ivan passed by within 50 km producing 50-mph 
winds, 5 tornadoes over the bay system with 70 mph gust winds and wave surf of 1.7 
meters (NOAA data). This extreme example is not lonesome and monthly wind speed 
increase between June and September is mainly due to tropical storm systems hitting that 
region. 

A study based on 110 year record (1886-1996) of tropical storm entering the 330 
km threat radius around Panama City with specific interest on Hurricane Eloise in 1975 
and Hurricane Opal in 1995 conducted by the Navyi shows that no location on St. 
Andrew Bay at Panama City is suitable as a haven for U.S. Navy ships in a hurricane 
threat scenario. The low-lying terrain of coastal Florida offers no protection from 
hurricane effects and the bathymetry of the Gulf of Mexico coast adjacent to Panama City 
is favorable for storm surge generation (vertical rise in the motionless water level near the 
coast caused by reduction of atmospheric pressure and wind stresses on the water 
surface) whenever strong on-shore winds, such as might occur with a tropical cyclone 
passage, are present. Furthermore, the wave setup (superelevation of the water surface 
above storm surge level due to onshore mass transport of the water by wave action alone) 
is enhanced by exceptional high waves occurring in such a particular period. A further 
complication results from the relatively slow speed capability of the Navy vessels 
(usually belonging to Mine Warfare Fleet) visiting Panama City, and their intolerance for 

1 available online at https://www.cnmoc.navy.mil/nmosw/tr8203nc/panamac/text/frame.htm the 20- 
May-05 
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seas in excess of 4 to 5 meters. During the largest events that hit Panama City, all these 
phenomena summed up to reach 3.5 to 4 meters and worst case scenario predicts more 
than 6 m. Figure 10 (Fig. XXVIII-6 of the aforementionned study) presents how these 
events are distributed within the year. 
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Figure 10, Annual distribution of the 126 tropical storms and hurricanes pissing 
within 330km of Panama City between 1886 and 1996 
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IV. 


HYDRODYNAMIC MODEL 


A, OVERVIEWS 

The model used in this study is a three-dimensional, time dependent 
hydrodynamie model, the Water Quality Mapping and Analysis Paekage (WQMAP) 
developed by Applied Seienee Assoeiate, Ine. out of Narragansett, Rhode Island. It is 
already being used by many institutions as NAVO, oil eompanies or eoast guards 
(Ireland, Spain and USA). Continuity, momentum, salt and temperature transport 
equations are solved after applying a generalized orthogonal eoordinate transformation on 
the horizontal and a sigma eoordinate transformation on the vertical. These resultant 
equations are then decomposed into an exterior mode (barotropic or vertical averaged 
equations) solved implicitly and an interior mode (baroclinic or vertically structured 
equations), solved semi-implicitly. Some approximations are made in the model 
formulation such as the hydrostatic balance, and Boussinesq approximation. 

This model solves the three dimensional conservation of mass, momentum, and 
energy equations with boundary conforming grid system and can be applied to both 
estuarine and littoral regions. The eddy viscosities can be both specified by the user or by 
turbulence closure using the turbulent kinetic energy (TKE) model. Output of the TKE 
model can be used in conjunction with a prescribed mixing length to determine the 
vertical eddy viscosities. The model outputs timeframe fields of surface elevations and 
velocities. The model can be forced with freshwater inflow, surface elevations, and wind 
fields. 

In the exterior mode, the Helmholtz equation, given in terms of the sea surface 
elevation, is solved to ease the time step restrictions normally imposed by the 
aforementioned gravity wave propagation. In the interior mode the flow is predicted by 
an explicit finite difference method, except that the vertical diffusion term is treated 
implicitly (Madala and Piaczek, 1977). The time step generally remains the same for both 
exterior and interior modes (Spaulding et al, 1999). Hence, the user can specify the 
variable (currents, temperature, surface elevation and salinity) over a certain simulation 

period and then the results/ predictions are presented in either color contours or vectors. 
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B, COORDINATE SYSTEM 

The basic foundation of the approach is to construct grids such that all domain 
boundaries are coincident with coordinate lines. This can be done using tensors which 
permit the model to transform the boundary fitted grid to a numerical grid employed for 
spatial discretization utilized in an Arakawa C Grid (Chu et al, 2004). Even though the 
transformed set of governing equations is considerably more complex than the original 
set, the transformed boundary conditions are specified on straight lines and the coordinate 
spacing is uniform in the transformed space. Orthogonal and conformal curvilinear grids, 
as well as simple stretched rectangular grids, are special cases of the generalized 
boundary conforming grid system (Spaulding et al, 1999). 

The two metric tensors to transform orthogonal position ((p,0) into conformal grid 
position (^,v|/) are defined by the following equations: 


^11 = 




cos^ 0 + 




(1) 


S22 


d(p 

dy/ 


cos 9+ 


( 2 ) 


y 

The user defines every single boundary point and then WQMAP numerically 
solves the coupled, elliptic coordinate transformation equations to determine the interior 
grid points within the model domain. 


C. GOVERNING EQUATIONS 

The hydrodynamic model included in the system solves the three dimensional, 
conservation of water mass, momentum, salt and energy equations on a spherical, non 
orthogonal, boundary conforming grid system and is applicable for estuarine and coastal 
areas (Muin, 1996). The two dimensional equations can also be predicted by vertically 
averaging the above over the water column. 
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A sigma stretching system (number of layers constant whatever the depth is) is 
used to map the free surface and bottom to resolve bathymetric variations. The basic 
dynamical system is represented by Kantha and Clayson (2000) with the following 
governing equations. 


The momentum equations in the and p-directions are given by 


duD 


- + - 


1 


^jgngn 
_ gD 


diuvD^A 5(V^) 

' - + uvD 
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dvO) ^ 

+ —- fDu 
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where r\, a, f, R, po and Ay successively represent the surface elevation, the vertical 

Z — 77 

stretching coefficient (a = -—), the Coriolis parameter (f = 2f2sin0), the Earth radius, 

H + 7J 

the mean water density (1025 kg/m^) and the vertical eddy viscosity. The vertical 
averaging of the afore equation leads to 
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+—(Tj-i)+A,pV^U 

Po 


where Ah, x'^ and x’’ sueeessively represent the horizontal eddy viseosity, the wind stress 
and the bottom stress. 


The continuity equation is given by 



where ro represents the vertical velocity normal to the sigma level. Finally, the transport 
equation for any substance S is given by 


dS ^ u dS ^ V dS ^ dS _ 1 1 ^ ^ fz) 

dt^R^di/^da P^gX "d^y R^gX "dy/y If day ^da) 


where Dh and Dy stand for horizontal and vertical diffusivities. 

D, STABILITY CONDITIONS 

The salt and temperature transport equations are solved by a simple explicit 
technique, except for the vertical diffusion term that is solved by an implicit scheme to 
ease the time step restriction due to the small vertical length scale. The advection term is 
solved using an upwind scheme (first-order accurate) that introduces artificial diffusivity. 
The horizontal diffusion term is solved by a centered in space explicit technique. The 
time step limitation satisfying the stability criteria is given by: 
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1 


( 9 ) 


^t< 



E. BOUNDARY CONDITIONS 

The flow normal to the shoreline is set to zero and there is no flow through the 
water surface or the bottom where the momentum equations successively reduces to: 


A^( du dv^ 



( 10 ) 


f du dv^ 
D \d(j’ d<J^ 





( 11 ) 


where Cd, Pa and W respectively being the drag coefficient, the density of air and the 
friction velocity. 

There is no transport of salt at closed boundaries as the shoreline is supposed 
impermeable. At the open boundaries, the concentration is specified during the inflow, 
using the characteristic values for Gulf of Mexico (temperature was set at 21“ C and 
salinity at 35 psu). The same temperature was used at the intracoastal open boundaries 
and the salinity was set to 22 and 20 psu at West and East Bay respectively according to 
Blumberg and Kim (2000). 


F. CONFIGURATION 

The results from Blumberg and Kim (2000) being considered, in this particular 

study, we ran the model with a rectangular grid in its baroclinic mode. The number of 

layers was set to 10, every layer representing one tenth of the averaged water column 

depth within every cell. The model is sensitive to bathymetry changes occurring between 

two nearby cells. To avoid it blowing up, smoothing passes were applied on top of the 

grid averaging and each cell was assigned a minimum 2m-depth to avoid a large 

computational flux increase. This process led to significant depth changes over the whole 

bay. As the model conserves the total mass, the resulting speeds computed in the different 
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cases correspond to the smoothed depths. Careful attention must be taken if these values 
have to be used for operational applications. 

Since the temperature is well mixed through the water column, we considered it 
as being constant in order to ease the computational effort. The 3-D choice is obvious and 
will not be discussed; the choice of a rectangular grid was driven by computational time 
restriction. For satisfying the CFL condition, the time step had to be set to 0.1 minute, 
and hence each run was requesting 17 hours. 
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V. 


DATA 


A, RIVERS 

Overall, there are 15 sourees of freshwater that are used as part of the model 
foreing. They are identified in Fig. 11 and tabulated in Table 3. The major souree of 
freshwater to the bay is Deer Point Lake. Daily-based flow rates were ealeulated using 
measured weir heights at the dam itself. The long-term average flow is 27m /s (Musgrove 
et ah, 1965). Freshwater also enters the bay at many other loeations. Unfortunately, these 
sourees are largely ungauged even through the flows are often substantial, espeeially 
during times of intense preeipitation. Table 2 represents the flow rate estimation 
eomputed by Blumberg (2000). 



sub-basin 

■■ 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

m^/s 

0.5 

0.8 

- 

1.3 

1.4 

0.4 

1.4 

1.7 

1 

1.2 

1.2 

1 

3.8 

- 

0.7 

0.9 

38.1 


Table 3. Calculated mean flow (Blumberg and Kim, 2000) 
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Nevertheless these data are not as erueial, and thereby as aeeurate as they would 
have been for mass transfer eomputation. The effeet of their fluetuation is, in this study, 
what is sought. 

B, BAY ENTRANCES 
1. Sea elevation 
a. Tides 

In this ease also, the aeeuraey was not a key request for running this 
model. Data were eolleeted at Panama City Beaeh where NOAA’s National Oeean 
Serviee maintains a station providing among others water elevation in m above MLLW 
(Mean Lower Low Water) every 6 (or 60) minutes. As we eonsidered the flow entered 
the bay through two passes, the evolution of each tidal constituent was examined as they 
propagate from East Pass onto West Pass on the continental shelf The distance 
separating the two channels is 11.5 km. As a result, the phase shifts were of the order of 
the degree (1.1° for diurnal constituents and 2.5° for semi-diurnal constituents) and 
corresponded to a 5 minute-time lag. Moreover, the slight difference on the bathymetry in 
front of both entrances did not involve a large magnitude difference (0.1% for semi¬ 
diurnal constituents - K1 or 01- up to 2% for shallow water constituent - M4 - which 
only contributes little (1.1 cm). At the end, these modifications were ignored. At last, even 
if the time series used at both locations were issued from Panama City Beach, 30km north 
of East Pass, we would not have had to correct them because the results from these series 
had not to be compared with any measured data set. 

However, NOAA/NOS does not maintain gauge at the three open 
boundary locations and therefore the elevation time series were got using Wxtide32 
software2 (also called Tides and Currents) which is a free tool providing pure tidal 
elevation all around the world. The following table represents from left to right the 
Eourier decomposition of elevation gauge time series collected at NOAA station on 
Panama City Beach, the tidal constituents used by NOAA for predicting tides at Panama 
City Beach and the Eourier decomposition of tide predictions from Wxtide32 software. 
Numbers in parenthesis represents the 95% confidence interval for the computed values. 

2 Available online at www.wxtide32.eom the 20-May-05 
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Panama City Beach 

collected data study 

NOAA prediction data 

Tide and current software 

amplitude [cm] 

phase [deg] 

amplitude 

phase 

amplitude 

phase 

K1 

15.89 (0.3) 

296.2 (1) 

14.5 

286.7 

15.53 (0.1) 

313.9(1) 

01 

15.84 (0.3) 

284.8(1) 

14.1 

284.5 

15.95 (0.1) 

302.3(1) 

Q1 

3.4 (0.3) 

270.7 (5.5) 

3.1 

273.4 

3.5 (0.1) 

295 (4) 

M2 

3.4 (0.1) 

287.4 (1.5) 

3.4 

277.1 

2.3 (0.1) 

329.5(1) 

S2 

2 (0.1) 

303.1 (2.5) 

2 

274.5 

0.7 (0.1) 

325.2 (3) 

form ratio 

5.90 

— 

5.30 

— 

10.50 

— 


Table 4. Main tidal constituents comparison 

The major difference between the data from Tides and Current software 
and NOAA data lies in the semi-diurnal constituents which seem to be underweighted. 
However the resulting error remains low, being considered the low amplitudes of both 
M2 and S2. A comparison of both NOAA gauge and Wxtide32 software time series 
presented on the subsequent figure shows that differences are larger during spring tides. 



Days in 2004 


Figure 12, Time series comparison between NOAA and Wxtide32 data set 
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Finally, as NOAA could not provide time series for the lateral open 
boundaries, they were foreed with Tides and Currents data. To be eonsistent, we also 
decided to force the Gulf open boundaries with time series provided by this software 

b. Residuals 

The residual elevation composed from storm surge, wind-induced waves 
and wave setup was obtained from the NOAA gauge measurements after filtering the 
tidal elevation out. The times series depicts a rather uniform shape, except during the 
hurricane event in September. 



The residuals provide an average 7.7 em setup, its standard deviation 
being 13.6 cm. It is to notice that the hurricane event does not impact much the statisties. 
When disregarding it, the mean drops to 7.3 cm and the standard deviation to 11.8 cm. 
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These residuals were put on top of the tidal elevation when foreing the last run Gulf open 
boundaries. 

2, Wind 

The wind files were also eolleeted at Panama City Beaeh and have been applied 
unmodified and uniformly over the whole domain. The station provides data hourly in 
standard weather eonditions and inereases the sampling to every 6 minutes when extreme 
weather does oeeur. WQMAP automatieally manages the sampling differenees. 

3, Temperature 

Despite the seasonal fluetuation of temperature, we ran the different eonfiguration 
setting a eonstant 2rC temperature value. As the temperature is uniform over the water 
eolumn, its infiuenee on density distribution is negligible. Foreing the fresh water 
temperature to the same value has probably introdueed the larger error. 

4, Salinity 

The different boundaries -West Bay, Gulf entranees and East Bay- were given 
salinity values -22, 35, 20 psu respeetively- in agreement with measurements eolleeted by 
Blumberg and Kim (2000). Finally, the rivers eells were foreed with pure fresh water 
whieh remains questionable as water enters the system through seepage. 
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VI. 


NUMERICAL EXPERIMENTS WITH WQMAP 


Different eases were featured with separate forcing to examine the special impacts 
caused by each of them. Three of them are described in this chapter occulting the run 
where the open boundaries were forced with residual elevation. It is enough to say that 
the results practically matched the ones observed in the case where open boundaries were 
forced by wind only. This being set, the following cases describe no forcing, tidal forcing 
and wind forcing situations. 

A. RUN 1: NO FORCING 

1. Initial conditions 

This case is featured with no wind and tidal forcing. The salinity is taken its 
climatological values at the open boundaries of the bay. The average values of runoff are 
assigned at all the river entries with pure fresh water. The salinity initial conditions for 
the basins were not crucial and several runs made from different values (32, 25, 15) 
ended up with the same final equilibrium state. The temperature was set to 2rC and 
surface elevation is automatically given a 0-value. 

2, Numerical results 

a. Salinity profiles 

The model is integrated from 01/01/2004 with no wind and tidal forcing. 
After one month, the model reaches its steady state with different zonal features. Fed by 
Deer Point Lake Dam, North Bay rapidly mixes and becomes fresh until it connects West 
Bay. St Andrew Bay presents a uniform 35 psu-salinity profile due to its connection to 
Gulf of Mexico. The salinity profiles of both East and West Bay are driven by the 
bathymetry contour. Being very shallow East Bay does not permit a lot of saline water to 
enter from St Andrew Bay and therefore, we end up with a rather fresh well-mixed 
profde of 4.8 psu. On the other hand, because West Bay is much wider and deeper, most 
of the salty water entering the system finishes its course westward and mixes the fresh 
input from Deer Point Lake Dam. The resulting feature depicts a uniform salinity profile 
of 17.5 psu. 
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A cross section view between West and East Bay points (Fig. 14) shows 
how salinity is distributed over the water column (Fig. 15). Two spots gather the 
conditions for a baroclinic circulation to take place and locate at marks A and B. They 
obviously stand into a transition zone between the salty inflow from the Gulf and the 
fresh runoff occurring farer inward. The resulting profiles (Fig. 16) portray a large 
stratification process. 
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Figure 16, Salinity profile at locations A and B (see Fig, 2) 


b. Velocity profiles 

The speed profiles only depend upon the salinity distribution as this run 
has no other foreing. Therefore, the speed is uniform over the water eolumn at all open 
boundaries but presents large change at locations A and B according with the salinity 
profiles described above. The velocity fields at surface and bottom (Fig. 17) clearly 
picture the baroclinic circulation at locations A and B. The velocity profiles at both 
locations are plotted in Fig. 18. 
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Figure 17, Velocity at (a) surface and (b) bottom (30° 09’ N - 85° 42’ W) 


averaged speed averaged direction 
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speed [m/s] direction [deg] 


Figure 18, Vertical variation of horizontal velocity: (a) speed at location A, (b) 
direction at location A, (c) speed at location B, and (d) direction at location B, (see 

Fig.2) 

c. Fluxes 

The computation of the flows at each open boundary describes how the 
water is flowing over the whole system. The total inflow represents 879 m /s (622, 202 
m /s at West and East Pass respectively and 55 m /s river runoff) roughly balanced by the 
water out-flowing at West and East Bay (538 and 334 m /s respectively). 

Accordingly with Eigure 18, the flow at locations A and B is two-layered 
and the system can be summarized as in Eigure 19. At locations A and B, dashed arrows 
represent the flow in the lower layer and the plain ones the flow in the upper layer. Notice 
that, despite the shear occurring at location A, the water always flows northwestward. 
The upper arrow symbolizes the fresh water input being aware that no location is 
associated with this arrow. It helps, however, the system to be in balance. All values are 
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in m /s and are not associated to standard deviation as the system is in steady state at the 
end of the run. By running the model in several river salinity values, we observed that 
increasing the river salinity by Ipsu reduced the inflow by 15 m /s. 



B, RUN 2: TIDAL FORCING 

1. Initial conditions 

This case is featured with purely tidal forcing at all open boundaries. They were 
initialized with different sets from Tides and Currents software. Both salinity and 
temperature were given similar values from the ones used in the first case. 

2. Numerical results 

a. Salinity profiles 

By driving the system with tidal forcing, we induced large fluctuations in 
speed and salinity profiles. The general features described in run 1 are still valid but the 
transition parts have been extended dramatically. On these spatial variations overlay 
temporal ones which are of the nature of tides. 
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First of all, the boundaries do not behave identieally. Beeause East Pass is 
too shallow, most of the water outgoing the system flows through West Pass. As a result, 
the salinity remains constant at 35 psu. West Pass, on the other hand, can see its salinity 
varying from 35 down to 21 psu, the larger fluctuation occurring in spring tides. 
Moreover, the water column remains well homogenized and a slight difference is 
perceptible during spring tides ebb (1 psu between surface and bottom). 

The two lateral boundaries portray two very different situations. East Bay, 
shallow and narrow, can hardly oppose the natural outward fresh flow. Hence, salinity is 
entering the system during small episodes with spring tides flows (Eig.20). 



Figure 20, Salinity time series at East boundary 


West Bay entrance exhibits much larger exchanges between the 
Intracoastal Waterway and the bay. In that case again, the water column appears to be 
well-mixed which emphasized run 1 observations that described the flow at boundaries to 
be barotropic. We can also notice how salt gradually diffuses inward during neap tides 
(less than 20 cm elevation fluctuations) (Eig.21). 
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Figure 21, Salinity time series at West boundary 


The following plots utterly illustrate the tidal pumping process due to tidal 
forcing. These two shots have been captured from two successive low and high tides 
during spring tides (Fig.22). The tidal pumping also results in large salinity changes in 
West Bay. During spring tides, the tidal pumping makes the water flow in the same 
direction over the water column even at locations A and B where the flow was initially 
baroclinic. So, the constant supply of West Bay with salt from St Andrew Bay stops 
during ebb flows and the water freshens, the salinity reaching its stable value during the 
subsequent neap tide. 


38 




















































































Figure 22, Salinity cross-section view during (a) ebb and (b) flood period 


The conditions observed at both West Bay and St Andrew open 
boundaries perfectly correlate these remarks (Fig.23). The salinity remains constant 
throughout the water column and varies from 35 to 32 psu at St Andrew opening and 
from 22 to 7 at West Bay boundary. Both panels represent salinity profiles at low and 
high tides measured at St Andrew boundary which explains why the salinity increases 
during ebb tide at West Bay boundary. 
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Figure 23, Salinity profiles at (a) West and (b) Gulf entrances 
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b. Velocity profiles 

The flow at the Gulf boundary controls the flow at the lateral boundaries. 
Nevertheless, the inward flow in never very important because of its natural tendency to 
go out of the system. The same observation prevails for the Gulf entrance with an inverse 
tendency. Finally, the mean flow over a year-long run almost equals the flow during 
flood periods (Fig.24). At these locations the standard deviation of the flow is almost one 
order of magnitude greater than the flow itself. The last feature to notice is that the 
vertical speed shear observed at Gulf entrance (Fig.24c) is not due to a vertical 
stratification but is likely to result from bottom interactions. 
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Figure 24, Vertical variation of horizontal velocity: (a) speed at West open 
boundary, (b) direction at West open boundary, (c) speed at West Pass, and (d) 

direction at West Pass, 


Hence, the general conditions described during run 1 at both locations A 
and B are no more valid. The subsequent plots (Fig. 25) describe how extreme can be the 
changes over a tidal period. On these schemes, the tidal pumping effects are obvious at 
both locations. Location A exhibits a particular feature in the sense that the mean surface 
salinity is greater than its value measured during high and low tides. 
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Salinity profiles at locations A and B (see Fig, 2) 
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Figure 25, 














Notice that the flow is no more baroclinie during ebb flows at location A 
(Fig.26 b). The standard deviation of the flux ean be tenfold greater than the flow itself, 
which is mainly visible at loeation A. In that case again, most of the salty water 
originating at West Pass flows towards West Bay which is much deeper and so make 
location A conditions fluctuate over a wider range. At location B, the tidal pumping 
influence results in a 90° shift of the flow in both the upper and lower layers but the shear 
remains (Fig.26d). Its lighter influence finally outcomes in a standard deviation of the 
order of magnitude of the flow itself 

averaged speed averaged direction 
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Figure 26, Vertical variation of horizontal velocity: (a) speed at location A, (b) 
direction at location A, (c) speed at location B, and (d) direction at location B (see 

Fig. 2) 

The subsequent plots (Fig.27) show how the flow propagates between St 
Andrew Bay and West Bay. The starting profile of the left panel stands in the solid black 
line of its right. The tide is high and the water column starts to mix up before the global 
salinity decreases with the ebb flow. Starting from the dashed blue line on the left side of 
the right panel, the salinity increases significantly in the lower layers and the water 
column stratifies all through the flood tide (the last black line pictures the beginning of 
the next ebb flow). These two plots show why the surface salinity during both high and 
low tides is less than its average. During the ebb flow, the surface salinity first raises until 
the water column has mixed and then drops. 
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Figure 27, Salinity profile evolution over a tidal period 
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As the flow induced by the horizontal salinity gradient is weak, the tidal 
impact drastically increases its fluctuation. Figure 28 presents speed magnitude and 
direction at locations A and B (mean and mean plus/minus standard deviation). 
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Figure 28, Vertical variation of horizontal velocity: (a) speed at location A, (b) 
direction at location A, (c) speed at location B, and (d) direction at location B, (see 

Fig.2) 


c. Fluxes 

In terms of fluxes, we obtained very similar averaged values to the ones 
computed in the run 1; this result was expected as tides should cancel out over a large 
period of time. If the resultant flow is still directed towards the same direction at all 
locations, it may vary significantly with tides. As a result, the standard deviation of the 
flow can be as large as 6 times the flow itself (at Gulf entrance). 
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C. RUN 3: WIND FORCING 

1. Initial conditions 

This case is featured with the wind forcing. Hourly wind data are collected at 
Panama City Beach from the NOAA gauge. The wind field has been applied over the 
whole area. This run should provide specific features due to wind forcing. Moreover, we 
paid specific interest over the hurricane period (15-17 Sep.) when winds were blowing at 
25 m/s. The open boundaries were forced by salinity only with no elevation. 

2, Results 

As the resulting wind obtained from NOAA reeords is rather weak and oriented 
outward (averaged speed and direction are 0.12m/s from NW), the global situation has 
not dramatically changed from the steady situation described in run 1. However, temporal 
variability can be large in the occurrence of severe weather systems or singular wind 
regimes. 

a. West Bay boundary 

When westerlies are blowing, the foree faces the natural outward flow at 
West open boundary which reverses when the speed exceeds 6 m/s. As a result, the whole 
flux is directed towards East Bay 


49 











b. East Bay and Gulf of Mexico boundaries 

Surprisingly, both East Bay and Gulf entrances open boundaries behave 
simultaneously and have opposite responses to wind forcing. Despite the orientation of 
the bay, the flow at East Bay boundary does not reverse when easterlies are blowing 
because they feed the system through East Pass. But when winds are blowing from north¬ 
east at a speed above 8 m/s, the flow overturns in St Andrew channel acting as a tidal 
pumping process. However, this was not detectable during the days forced by Hurricane 
Ivan winds. When over 20 m/s winds from south-east forced the system, the entire flow 
directed towards West Bay. 

c. Hurricane Ivan 

We focused on the period when Hurricane Ivan hit the north Elorida 
Panhandle. We considered this period to begin when winds reached a steady 5 m/s speed 
on the 15th of September at 8 PM and to end 2 days later after wind has felt under 5m/s. 
Meanwhile, the average speed was 15 m/s and its peak achieved 25 m/s in the morning of 
September the 16th during 4 consecutive hours (Pig.30). Even though strong winds all 
come from the same direction, a rather slight change can largely impact the whole 
circulation in the bay. 
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If we noticed that East Bay reacted to wind forcing as it would have done 
to tidal pumping, West Bay response to strong wind forcing largely depends upon its 
direction. Fig.30 - 31 depict two dramatically different situations as the wind shifted from 
105 to 150 between 05;00PM and 01:20AM. Oriented in the axis of East Pass, the wind 
drains the water westward. Because the current is so intense, the water enters North Bay 
avoiding the fresh water to flow from Deer Point dam. As noticed above, the water also 
enters the system through East Bay accordingly with the flow at East Pass (Fig.31). 
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Figure 31, Surface current the 09/16 at 05:00 AM 


Later on after the wind had shifted, the eireulation pattern does not 
correlate anymore the previous one. In Fig. 32, the wind is blowing from 150°. Because 
the circulation reacts with West Pass flow increase, the water now exits the system 
through East Bay. The water, pilled up by the winds, still enters North Bay and hence 
there is a large deficit of water inside West Bay balanced by the western open boundary. 
A twisting motion then takes place for the surface current to balance the wind stress. 
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Figure 32, Surface current the 09/16 at 01:20 PM 


The wind also impacted the system by increasing its salinity. The 
elevation in North Bay increased up to 50cm from a standard value of 6cm, containing 
most of the fresh water supply in the very northern end of the system. As Hurricane Ivan 
passed by, we could detect lumps of salty water alternatively diffusing into West Bay and 
East Bay from all three open boundaries. 
d. Fluxes 

As the mean wind blew from NW, the resulting impact is weak and 
decreased the fluxes a little. The standard deviation is not large because strong winds 
were always blowing in the same direction and hence the flow never really reversed. 
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VII. 


FORCING MECHANISMS FOR THE HYDRODYNAMIC 

SYSTEM 


The control run, a case featured with all forcing mechanisms described in the 
previous chapter, is used as reference. We then disturbed the forcing conditions to assess 
their influence on the circulation and salinity distribution. The perturbations we applied 
have been chosen such as they are reasonably realistic. The results from the beginning of 
April and the end of September are compared applying error measurements between the 
control run and every single disturbed run. 

A, STATISTICAL ANALYSIS 

Forcing mechanisms are detected using three statistics applied on y/ representing 
(m, V, S) for both the surface and bottom layers: relative mean difference (RMD), root 
mean square difference (RMSD), and relative root mean square difference (RRMSD). 
Let {Wref^¥sen) ^hc modcl outputs [i.e., (u, v, S)] for the reference and sensitivity runs. 


The three statistics are defined by 


RMD{x,y) 



RMSD{x,y) 




RRMSDix, y) = ■ 

V n 


( 12 ) 


B. CONTROL RUN 

This case is featured with all forcing at all open boundaries. The model is 
initialized from zero velocity, 25 psu salinity and 21°C temperature throughout the whole 
domain. The model is integrated with the previously described elevation time series sets 
and the wind field. The boundaries are given the standard salinity and temperature values 
used in the preceding cases. The river cells took their averaged runoff values. 

The next figure gives a synthetic vision of the fluxes at open boundaries and 
locations A and B. The fresh water input has been represented with a single arrow which 
does not match the reality but was necessary for the mass conservation balance. All the 
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values are in m^/s. The flux at loeation A and B has been split into an upper flow (plain 
line) and lower one (dashed line). 



C. FRESH WATER INFLUENCE 

In this ease, the river eells are featured with a doubled total flux. All the 
remaining eonditions remain unehanged. The inerease (100%) ean appear very large but 
only represent a 55 m /s extra input. 

The extra fresh water input slightly impaeted the dynamies in St Andrew Bay 
system. As a matter of faet, this 100% inerease only represents 8% of the total amount of 
water flowing inward the bay and as 20m /s out of the initial 55m /s of fresh water input 
originate in North Bay, the impaet in terms of veloeity or salinity are negligible at the 
Gulf entranees (RMSD values of 3.10' m/s for both u and v eomponents and RMSD of 
0.16psu for the salinity) and East Bay (RMSD of 2.10' m/s for u- and v-eomponents and 
RMSD of O.lVpsu for salinity). The values of the RRMSD average 5% elearly deseribing 
a weak impaet. It is elear from Fig. 35 and 36, that the water eolumn remains well mixed 
and that the influenee is mueh greater at West Bay open boundary than at St Andrew 
Bay’s. At West Pass, the salinity disturbanee only oeeurs during spring ebb tides whereas 
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the influence at West Bay exit can be sensed during the whole period. The impact on East 
Bay was insignificant and will not be further discussed. 

As the initial disturbance constantly applies by increasing the inflow, the resulting 
errors for u, v and salinity always have the same sign for a given location. 
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Figure 35, Time series of (u, v, S) differences between the control and river 
disturbed runs valid for both surface and bottom layers at West Pass. 
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Figure 36. Time series of (u, v, S) differences between the control and river 
disturbed runs valid for both surface and bottom layers at West Bay open 

boundary. 


1 

The values, all in m /s, represent the flux eomputed at eaeh time series loeation 
and the differenee between this run and the eontrol run is indieated in parenthesis. The 
differenees observed at open boundaries eorrelate with the extra fresh water input. The 
highest values loeate between West Bay and St Andrew Bay where most of the water 
flows and where the greatest fresh water sourees stand. 
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D. WIND IMPACT 

This case was featured without wind foreing, therefore the system was forced by 
river inflow and surface elevation. During the studied period, the averaged wind was 
blowing from 155 at a 4m/s speed. 

This run shows clear features happening at all open boundaries. It is clear that the 
lack of wind has few not to say any impact on speed in the bottom layer even at East and 
West Bay open boundaries where depths average 2m. The bottom layer salinity is only 
influeneed at west Bay for whieh the presenee of high fresh water supply in North Bay is 
likely responsible. In this ease, we can notice that the impact can either be positive or 
negative depending on the direetion the wind was blowing from. The RMSD values for 
the speed are of the order of 10' m/s and the RRMSD ranges around 15%. The influenee 
on salinity however is signifieantly different at lateral boundaries and at West Pass. West 
Pass is eharaeterized by a 0.24psu-RMSD corresponding to a 0.7%-RRMSD whereas 
West Bay open boundary time series show a 2.4psu-RMSD and a 20%-RRMSD. The 
RRMSD is even higher at East Bay and reaehes 30% with a 0.7psu-RMSD. The lower 
salinity observed in East Bay enhances the weight of the relative error. 
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In the surface layer, the impact is obvious from East Bay to West Bay. The lack of 
wind enlightens both tidal and residual signals which can be easily identified. The speed- 
errors induced by the absence of wind are often greater than the speed itself and are 
tenfold higher than the errors computed in the bottom layer. The RMSD of the speed u- 
component fluctuates around 0.25m/s at West Pass and West Bay which correspond to 
RRMSD being 460 and -150% respectively (the negative sign portrays an outflow). At 
East Bay open boundary, RMSD and RRMSD values are respectively 0.08m/s and 117%. 
These errors do not weight identically when scouting the salinity. At West Pass, the 
RMSD has a value of 2.8psu which stands for a RRMSD of 8%; at West Bay open 
boundary, RMSD and RRMSD have values of 7.6psu and 60%; at East Bay open 
boundary, these values are 5.2psu and 200%. At last, one must notice that the fluctuations 
can achieve values as high as 2m/s for the speed and 25psu for the salinity depending on 
the location. 
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Figure 38, Time series of (u, v, S) differences between the control and no wind 
forcing runs for the surface layer at East Bay open boundary. 
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Figure 39, Time series of (u, v, S) differences between the control and no wind 
forcing runs for the bottom layer at East Bay open boundary. 



Figure 40, Time series of (u, v, S) differences between the control and no wind 
forcing runs for the surface layer at West Pass, 

61 










.2 I_■_■_■_—u_9_■_ 

120 140 160 180 200 220 240 

time [Julian day] 


Figure 41, Time series of (u, v, S) differences between the control and no wind 
forcing runs for the bottom layer at West Pass, 



Figure 42, Time series of (u, v, S) differences between the control and no wind 
forcing runs for the surface layer at West Bay open boundary, 
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Figure 43, Time series of (u, v, S) differences between the control and no wind 
forcing runs for the bottom layer at West Bay open boundary. 


Figure 44 plainly depicts the surface influence of the wind. It is interesting to 
notice that despite the lower outward flux occurring in East Bay, the absolute difference 
equals the one observed in West Bay, clearly pointing the shallow waters wind driven 
circulation. 
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E. WIND FLUCTUATION 

This case was features with a modified wind time series. We applied on each 
component of the original wind time series a normally distributed random perturbation 
whieh mean is 0 and standard deviation 1 m/s. All the remaining conditions stayed 
unchanged when compared with the control run. 

Because, in that ease, the disturbanee was randomly applied, the eomparison 
between the control run and this one gives results far different from the ones obtained 
when the wind had been removed. First of all, the surfaee RMSD values are very weak 
for both speed and salinity being 0.6cm/s and 0.04psu at West Pass, 0.4cm/s and 0.14psu 
at East Bay and lem/s and 0.4psu at West Bay. These values are approximately 25 times 
lower than the ones observed in the previous comparison. The differenee between the 
surfaee and the bottom layers is, however, insignifieant as shown by Fig. 45 to 50 and 
RMSD values remain praetieally unehanged in the deeper layer although slightly 
reduced. 
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Figure 47. Time series of (u, v, S) differences between the control run and the 
wind disturbed forced run for the surface layer at West Pass, 
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Figure 48, Time series of (u, v, S) differences between the control run and the 
wind disturbed forced run for the bottom layer at West Pass, 
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Figure 49, Time series of (u, v, S) differences between the control run and the 
wind disturbed forced run for the surface layer at West Bay open boundary. 



Figure 50, Time series of (u, v, S) differences between the control run and the 
wind disturbed forced run for the bottom layer at West Bay open boundary, 
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The lack of difference throughout the water column probably stands in the fact 
that the minor increase or decrease applied on the wind may not have induce a significant 
wind direction shift which impact is predominant as observed in the previous chapter. As 
a matter of fact, the highest differences observed in the previous plots occurred when the 
wind speed had its lowest values. The computation of the fluxes did not depict great 
changes as the perturbation respected a mean of 0 and Fig. 51 clearly shows that point. 



F, TIDAL RESIDUAL 

This case was featured in order to evaluate the impact of the tidal residuals. We 
then forced the model with tidal constituents only, getting rid of the residuals, and with 
the original wind data set, all remaining conditions staying unchanged. During this 
period, the residual averages 7.34 cm. 

The first observation from Fig. 52, 53 and 54 is that the lack of residual impacts 
both the surface and the bottom layers the same way and the difference between both 
layers is slightly discernible. For u and v, the impact the absence of residual created is 

more evident at West Pass (Fig.53). We can observe that the remaining signal has not a 0- 

68 












mean value. Because the average residual was positive, its absence diminished the inward 
flux and hence the difference between the control run and the new one is positive. The 
RMSD at that location is 0.07m/s for an RRMSD of 126%. 
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Figure 52, Time series of (u, v, S) differences between the control and no residual 
forcing runs valid for both surface and bottom layers at East Bay open boundary. 
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Figure 53, Time series of (u, v, S) differences between the control and no 
residual forcing runs valid for both surface and bottom layers at West Pass. 



Figure 54, Time series of (u, v, S) differences between the control and no 
residual forcing runs valid for both surface and bottom layers at West Bay open 

boundary. 





The error indueed on the eomputed fluxes is huge. Beeause the residuals are 
fundamentally positive (wave surge, quasi permanent high pressure systems above this 
area), the pressure gradient has tremendously diminished when they have been removed. 
Notiee that the impaet is greater towards West Bay. 



G. TIDAL CONSTITUENT FLUCTUATION 

This ease was featured in order to analyze the aeeuraey weight of the tidal 
eonstituents of the time series. This ease gathers three runs in whieh the Gulf open 
boundaries were sueeessively foreed with modified elevation time series. In the first ease, 
we deereased the diurnal K1 amplitude from 15.53 down to 14.5om and in the seeond one 
we inereased the semi-diurnal M2 amplitude from 2.3 up to 3.4om. These two new values 
eorrespond to the ones used by NOAA for predieting tides in Panama City Beaeh. The 
third run used the gauge elevation eolleeted at the NOAA Panama City Beaeh station. All 
the remaining parameters remained unehanged when eompared with the eontrol run. 

The first two oases did not produoe great ohanges on the system. The small 

variation applied sueeessively on the diurnal K1 and semi-diurnal M2 tidal eonstituents 

did not really impaet the pioture obtained with the initial oonditions in the eontrol run. As 

no shift was applied on the oonstituent phases, the perturbation obviously oorrelates the 

71 












tidal elevation. The speed RMSD values remain weak (1.2em/s at West Pass, 0.5cm/s at 
East Bay and 1.5em/s at West Bay) and so do the salinity RMSD values (O.lpsu at West 
Pass, 0.15psu at East Bay and Ipsu at West Bay). It is interesting to notice that, even if 
the disturbances were applied at Gulf entrance open boundaries only, the greatest errors 
were observed at West Bay open boundary. It emphasizes the fact that most of the inward 
flux flows towards West Bay where the shallowness enhances the speed shift. The speed- 
RMSD are constant at both West and East Bays through the water column and present at 
West Pass a 40% difference when comparing the surface layer to the bottom layer, 
salinity RMSD remaining unchanged. The RRMSD values however represent 20% in 
both layers, which tends to prove that the difference between the two layers is likely to be 
related to the bottom friction. 



Figure 56. Time series of (u, v, S) differences between the control run and the 
tidal constituent disturbed run valid for both surface and bottom layers at East Bay 

open boundary. 
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Figure 57. Time series of (u, v, S) differences between the control run and the 
tidal constituent disturbed run valid for surface and bottom layers at West Pass, 
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Figure 58. Time series of (u, v, S) differences between the control run and the 
tidal constituent disturbed run valid for both surface and bottom layers at West Bay 

open boundary. 


-5 

Over the five month period, the flux fluctuation is everywhere less than 3 m /s. 
The next figure, representing fluxes, is by the way valid for both changes applied on K1 
and M2 tidal constituents. 
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VIII. 


HYDROCHEMICAL PROCESSES AND CHEMICAL THREAT 


Besides WQMAP, ASA also built model applieations used as rapid response 
models named OILMAP, SARMAP or CHEMMAP after Oil, Seareh and Reseue or 
Chemieal Mapping and Analysis Program sueeessively. We ran the three dimensional 
version of WQMAP whieh prompted us to use CHEMMAP in order to evaluate the 
baroelinie influenee on aeeidental release. Eaeh ease will be featured with ehemieal 
eompound being released at surfaee and near the bottom to estimate the effeets on 
dispersion. We will also eonduet the release at different period of time to assess the 
influenee of foreing eonditions on the dispersing proeess. 

A, CHEMMAP 

CHEMMAP is a ehemieal diseharge model designed to prediet the trajeetory, 
fate, impaets and biologieal effeets of a wide variety of ehemieal substanees three- 
dimensionally in both marine and fresh water. The CHEMMAP model system is made up 
of several integrated eomponents among whieh the ehemieal database, the ehemieal fate 
model are of outmost importanee. 

1, The chemical database 

The database eneloses all the physieal/ehemieal properties required for simulating 
transport and fate of the spilled material. CHEMAP uses either the Chemieal Abstraet 
System (CAS) registry number or the ETN number to referenee eaeh ohemieal. The 
ehemieal state ean be defined as solid, liquid or gas, dissolved or not in aqueous solution 
or in solvent. As ehemieal properties do vary with temperature, there all referred to their 
value at 25°C. 

2, The chemical fate model 

The spill model itself prediets the movement of ehemieals on the water surfaee 
and the distribution of ehemieals in the environment in the water eolumn, in sediments 
and on the shoreline. The model simulates spreading, dispersion, evaporation, 
entrainment, dissolution, sedimentation and degradation proeesses. The model is 
initialized with the user-speeified time, loeation and duration release of seleeted material. 
Its state and solubility are driven the initialization algorithm. In the ease of a highly 
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soluble chemical, the chemical mass is initialized in the water column in the dissolved 
state. All the different cases are defined in McCay and Isaji (2004). 


The mass is transported in three-dimentional space and time, relying on 
environmental data such as wind drift and currents. The vertical velocity of chemical 
components is computed with Stoke’s law, which relates the terminal settling or rise 
velocity of a smooth, rigid sphere in a viscous fluid of known density and viscosity to the 
diameter of the sphere when subjected to a known force field (gravity). The equation is: 


w- 


2 Pwater^ ^chem S 


V. 


(13) 


where plater the densities (kg/m ) of the chemical component and the water; 

fchem is the radius (m) of the chemical component; is the water viscosity (Ns/m ). 

The currents in the surface water layer (surface drift currents) are calculated with 
hourly wind field. The surface drift currents have the magnitude ranging between 1 and 
6% of the wind speed and are in a direction 0-30° to the right (in northern hemisphere) of 
the down-wind direction. In this study, the horizontal turbulent diffusion coefficient 
ranges between 1 and 10m /s and its value was set to 1 m /s in the different cases. The 
vertical turbulent diffusion coefficient ranges from 0.00001 (for stratified ocean) to 0.001 
mVs (for well mixed ocean). It is set to be 0.00001 m^/s in this study. 

The chemical is modeled using a Lagrangian approach in which spillets are 
tracked in both space and time. All the phase transfer rates are computed at each time 
step. This type of model is frequently used to trace back harmful substances to their 
sources and constitutes a valuable tool in the identification of environmental polluters. 

3. Environmental data and stochastic model 

CHEMMAP also includes an embedded Geographic Information System (GIS). 
The GIS is used to store, display and analyze any type of geographically referenced data. 
GIS also includes data which are not necessarily used by the spill model, such as place 
names, critical habitats for fish and wildlife, spill response equipment, shipping lanes. 
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and real-time spill observations, but they are often helpful in analyzing and interpreting 
model results. 

The stoehastie model randomizes at a given loeation spill date and henee wind 
and current conditions. This model output consist in statistical maps defining, for 
example, the probability of exceeding a concentration threshold or the time when the 
concentration first exceeds the threshold. 

B. CHEMICAL COMPONENTS 

The choice of the chemical components used in the afterward scenarii ensued 
from our will to analyze the baroclinic impact on dilution and spill. Therefore, instead of 
focusing on the lethal power of the chosen component, we picked up components which 
tend not to evaporate promptly. Since this thesis is unclassified, we are not going to relate 
the results to any relevant threat scenario. However, our statements can be transposed to 
all similar chemical components which are in the field of interest in national security 
concerns. 

1. Physical information 

Ethylene glycol, OH-CH 2 -CH 2 -OH, is a non volatile, soluble sinker, its density 
being 1.1132 kg/m^. It is referenced as CAS 107-21-1 or UN 8027. Ethylene glycol 
is a clear, colorless, slightly syrupy liquid at room temperature. It is odorless but has a 
sweet taste. 

2. Industrial use 

It is used to make antifreeze and de-icing solutions for cars, airplanes, and boats; 
to make polyester compounds; and as solvent in the paint and plastics industries. 
Ethylene glycol is also an ingredient in photographic developing solutions, hydraulic 
brake fluids and in inks used in stamp pads, ballpoint pens, and print shops. 

3. Exposure effects 

Eating or drinking very large amounts of ethylene glycol can result in death, while 
large amounts can result in nausea, convulsions, slurred speech, disorientation, and heart 
and kidney problems. Eemale animals that ate large amounts of ethylene glycol had 
babies with birth defects, while male animals had reduced sperm counts. However, these 
effects were seen at very high levels and would not be expected in people exposed to 
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lower levels at hazardous waste sites. Ethylene glycol affects the body's chemistry by 
increasing the amount of acid, resulting in metabolic problems. Similar to ethylene 
glycol, propylene glycol increases the amount of acid in the body. However, larger 
amounts of propylene glycol are needed to cause this effect. 

Its primary hazard is the threat to the environment. Immediate steps should be 
taken to limit its spread to the environment. Since it is a liquid it can easily penetrate the 
soil and contaminate groundwater and nearby streams. 

4, Safety recommendations 

The EPA has set a drinking water guideline for ethylene glycol of 7,000 pg in a 
liter of water for an adult. The American Conference of Governmental Industrial 
Hygienists (ACGIH) recommends a maximum level of 127 milligrams of ethylene glycol 
per cubic meter of air (127 mg/m^) for a 15-mmute exposure. 

C. TRANSPORT PATTERN 

In all cases, 10 tons of the aforementioned chemical constituent are released 
within 10 hours with an original plume thickness of 0.5m. The choice of the date did not 
respond to any particular need. The release locations, however, were selected in order to 
maximize the tidal pumping effect and minimize the vertical mixing process. 

1. Release at (30° 08’ 45.5”N, 85° 40’ 46.8”W) on June 1*‘ at 12am, 

The chemicals are released at (30° 08’ 45.5”N, 85° 40’ 46.8”W) 12 am on June 
H*. We analyze differences between surface and bottom release. As ethylene glycol is not 
volatile, the evaporation process is weak. Because of the tidal pumping, the chemical 
dilutes far away from its source in both West Bay and East Bay directions (Eig.61). Even 
if this location is closer to East Bay, it essentially diffuses westward accordingly with the 
mean winds (from 155° at 4m/s) whereas flow tends to push the water eastward. The 
maximum dissolved concentration rarely exceeds Img/m beyond 85°35’ toward East 
Bay and occurs between 4 and 10 days after the release. The same peak reaches a 
position as far as 8.5km south west of West Pass on the inner shelf 4 days after the 
release. A maximum concentration of 32000 mg/m is obtained 1 hour after the release 
started (Pig.62). Pigure 63 shows that the ethylene glycol is immediately mixed into the 
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water column. Its decreasing rate is 0.3024 per day at 25°C and after 3 weeks ethylene 
glycol has disappeared from the bay. 



Figure 60, Wind field over the domain between the Ol-Jun-2004 12:00am and the 

21-Jun-2004 11:00pm 


81 



Concentration (mg/m3) 



Figure 61, Spillet dispersion after 3 weeks 
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Figure 62, Maximum concentration for dissolved ethylene glycol 
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Mass Balance for ethylene glycol 

Surface Water Column Ashore Evaporated 



Time (days) 

Figure 63. Mass balance for phenol release at location tp = 30° 08’ 45,5N -G = 85° 

40’ 46.8W 


It is difficult to correlate the signal between two different loeations as the plume 
goes baek and forth with tides and as the eoneentration depends upon the depth at eaeh 
loeation. As it propagates westward from the release point to the end of West Bay, the 
shape of the signal varies signilieantly as shown in Fig. 64, 65 and 66. We ean also notiee 
that, the shape smoothes with time (only the largest peak ean be traeked) but when the 
depth is shallow enough, we deteet again higher peaks whieh were present nearby the 
release loeation. It is elear from these plots that the water depth influenees the dissolved 
eoneentration between the last two loeations (the depth deereases from 10 to 2m). The 
eoneentration naturally deereases as we move away from the souree. 
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Figure 64, Signal in St Andrew Bay 


Maximum Dissolved Concentrations 


15.0 

12.5 
10.0 

7.5 
5.0 

2.5 

0.0 

Figure 65, Signal at location A 
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Maximum Dissolved Concentrations 

85.8409W, 30.2785N for ethylene glycol 
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Figure 66, Signal at West Bay open boundary 


Figure 67 enhances the effects a semi enclosed small creek can produce. With 
southeastern winds, the chemical piles up in the northern part of West Bay where it 
remains and accumulates until being diluted into West Bay by changing wind conditions. 
We can therefore find lumps of pollutant in very small areas which can more severely 
impact the environment. 
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Figure 67, Signal in the north of West Bay 
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2, Wind influence study 

If the tidal influence in the spreading of the plume can be viewed during the 
simulation, the impact of the wind is less obvious. Two cases are designed to show the 
wind effect: (1) reversing wind direction and (2) eliminating wind. When the simulation 
was done without wind, most of the pollution remained contained between points A and 
B defined in the previous chapters during the first two weeks (the 15* of June 
corresponds to the highest tidal range after the release). After that period, a large part of 
the chemical diffuses into West Bay. Figure 68 also shows the following features. As the 
pollutant transits within the mean flow in West Bay, it seems to pile up in different places 
in East Bay. This statement only characterizes that particular snapshot and does not 



Figure 68, Spillet dispersion after 3 weeks without wind 


When the wind field is reversed, it is not surprising that the pollutants mainly 
diffuse eastward (the average wind now blows from WNW at 4 m/s) and outward the bay 
system (Fig.69). 
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Figure 69, Spillet dispersion after 3 weeks with reversed wind 


3, Release time influence using the stochastic model 

The dependence on the wind field appears to be of importance. None of these 
cases, hence, is representative of the average impact a chemical release would have on 
the bay even though the release location remains unchanged. We then applied the 
stochastic model featured with 50 randomized dates of release taken between the of 
June and 31®‘ of August. Figure 70 represents the run number of the worst case scenario 
for each cell (each run is featured with different release date, then different wind and 
current data but same location and same amount of chemical). The worst case scenario is 
defined for giving the maximum value encountered at each cell during all the runs. Thus, 
the large discrepancies observed on figure 68 lights the importance of the chosen date of 
release. To each tone of grey correspond 5 consecutive runs (even if they are not related 
to each other). 
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Figure 70, Run numbers for worst case scenario 


Fig. 71 shows a severe ease with a maximum dissolved coneentration and two 
major features. First, the pollutant predominantly diffuses westward whieh eorrelates the 
averaged wind encountered during this period of the year. Second, the pollutant tends to 
concentrate in shallow motionless waters on the edges of West and East Bay. Note that 
the two patches of lower concentration in St Andrew Bay correspond to the deepest 
locations of the bay. 
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Figure 71, Worst case maximum dissolved concentration (mg/m^) with spill 

released in St Andrew Bay 


Therefore, defining the probability a threshold to be overshoot at each grid cell is 
important. This threshold represents a short term exposure limit (STEL) or an 
immediately dangerous for life or health (IDLH) limit. The model provides such a 
probability being given in percentage of runs during which this threshold has been 
exceeded. If we only consider the possibility of the bay system to be polluted (threshold 
is given a 0-value), then this percentage is 100% for the whole bay except in the northern 
part of North Bay and the eastern one of East Bay where probability never drops under 
50%. The knowledge of such a probability is vital for deciding the deployment of 
antipollution devices. But of even furthest importance is the time given before this value 
can be reached. Eig. 72 then represents the minimum time requested for exceeding the 
threshold at each grid cell. Hence, for that particular release location, all St Andrew Bay 
will be polluted within 12 hours, within 24 hours the pollution will extend from point A 
to the first third of East Bay. It will last over 3 days for the pollutant to reach West Bay 
open boundary and 10 days for East Bay’s. 
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Figure 72, Minimum time (hours) to exceed threshold 


4. Spill location influence study 

It is obvious that the time for the pollutants transporting outward the bay mainly 
relies on the winds and source location. If the application of the stochastic model 
somehow handles the wind factor, the release location must be shifted to other places in 
order to evaluate the relative weight of this factor. We then applied the same concept to 4 
different sites, point A, point B, the center of East Bay and the center of West bay which 
were supposed to depict different pollutant propagation features accordingly with the 
hydrodynamic model results. 

Because the flux originating from Gulf of Mexico predominantly flows westward, 
a release point located between St Andrew Bay and West Bay or even more inside West 
Bay is likely not to deeply impact East Bay. The pollution will only reach the end of East 
Bay after 15 days if the pollutant is spilled at point A, which only stands 8km away from 
the previous spot. As this chemical is not volatile, it does not evaporate and its mass is 
roughly conserved into the system until its natural decay acts. West Bay is much 
shallower than St Andrew Bay, thus the small concentration decrease close to the Gulf 
entrance involves a large pollutant mass transfer into West Bay. 


90 








■)20 

■)30 

■)I0 

■)50 


■)iW 

■)70 

■)S0 


■)30 


■) 1001 



Figure 73, 


Worst case maximum dissolved concentration (mg/m^) with spill 
released at point A 
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IX. 


CONCLUSIONS 


In this thesis, we first study the baroelinieity and foreing meehanism for a 
hydrodynamie model in St Andrew Bay, Florida. The un-nested trait of this model 
required the open boundaries to be loeated elose to the study domain in order to apply a 
uniform wind profile with some eonfidenee. If plaeing the boundaries more offshore 
should provide better results in term of eireulation at Gulf entranees essentially, it would 
also induee problems we did not faee. Buoys deployment would be required in order to 
get sea level data and all the open boundary eells eould not be foreed with the same data 
set beeause of the phase shift indueed by the stretehing of the boundary. The use of 
smoothing proeess in very narrow fast bathymetry-ehanging bays dramatieally impaets 
the absolute values of both speed and overall fluxes. Therefore, speed eurrent values, 
when required for operational purposes (diving operation for example), must be used 
with eaution. The assessment of this proeess, whieh was neeessary not to inerease the 
number of eells signifieantly and henee deerease an already very small time step, should 
request further study. 

Despite these shorteomings, the baroelinieity of the eireulation was pinpointed 
and the ehannels joining St Andrew Bay to both West and East Bays have elearly been 
identified as being two-layered. Notiee that, even though no rivers are present on seene, 
the freshwater input from underground seepage was eharaeterized in the model as 
originating from rivers. It is then reasonable to estimate that the faet the two-layered flow 
is eonfined in two areas takes its origin in this bias. This result alone is suffieient for 
using a 3D baroelinie model. However, the magnitude of the freshwater input does not 
really infiuenee the flux itself neither did the aeeuraey of the wind field or of the tidal 
eonstituent amplitude. Nevertheless, the impaets of wind on vertieal stratifieation and 
those of tidal residual on the flow are important, partieularly for operational purposes, if 
all foreing meehanisms are not taken into aeeount. 

The baroelinieity was not however an outmost feature in the CHEMMAP 
applieations. Beeause the ehemieal data embedded provided substanees whieh were 
either immediately dissolved into water after their release, like ethylene glyeol, or into 
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the sediment, like tetraethyl lead, or even too volatile for remaining into the water, we 
eould not study the differenees implied by releases at surfaee and bottom. We plainly 
proved the furthest importance of the wind in the pollution drift, particularly in St 
Andrew Bay system where the shallowness gives an outstanding weight to the wind in 
driving the flow. Because the water flows inward at Gulf entrances and splits 
asymmetrically mostly towards West Bay, the dependence of the chemical dispersion on 
the release point is strong and East Bay is hardly impacted unless the pollution takes 
source in the bay itself or when winds are blowing eastward. The chemical clearly 
propagates westward during the simulation with a concentration generally damping away 
from the source except in very shallow zones where some accumulation can occur with 
favorable winds. These accumulations are finally responsible of secondary pollution 
events. It is, at last, obvious that these results were biased, again, by the bathymetry 
smoothing process. However, it is reasonable to conclude that, as East Bay presents a 
rather uniform depth, this bias impact was not so crucial in this area and also that most of 
the pollution drift relied on the wind. 

Eurthermore, the description of different salt diffusion processes during ebb and 
flood tides besides the observation that spring tides were globally causing a general 
freshening of the system by pumping fresh water while neap tides let salt to diffuse from 
the saltier St Andrew Bay towards both East Bay and particularly West Bay constituted 
two central resulting outcomes. The tidal impact study finally described an estuarine 
circulation with imbalanced ebb and flood periods. 
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